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The problem of approximating a band-limited function, H(t), by a sum 
of cosines arises in the design of phased array antennas. Three methods of 
synthesis are presented for establishing such designs. Error formulae are 
deduced for each method, including a new error formula for Tchebycheff 
quadrature. The existence of grating lobes is proved, and lower bounds for 
their location are developed. 

I. INTRODUCTION 

This paper is concerned with the problem of approximating a function 
H(t) (-co < t < oo)bya cosine sum of the form 

1 N 
S w (t) = rr= £ cos txj , ^ xi < Xi < ■■■ < x N ^ 1. (1) 

JSI ;=i 

The synthesis of array antennas is an application of the problem of this 
paper. Let isotropic radiating elements of strength \/N be located 
along the .r-axis at the points Xj{\ ^ j ^ N) providing planar radiation 
of wavelength X, and let 6 designate the angle between the positive 
?/-axis and a line passing through the origin and a far-field point, then, 
setting 

t = 2 JL^l, (2) 

A 

the far-field radiation pattern of the linear array is given by SN(t). 
The requirement that all the coefficients of the sum in (1) be equal 
generally stems out of the use of identical radiating elements, and out 
of the desire to employ identical feed for each element. 

The function H(t) represents the desired far-field radiation pattern; 
it will be required to satisfy the condition 

H(t) = / F(x) cos tx dx, (3) 
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for some function Fix) 6 L(0,1), and the normalization condition 

H(0) = 1. (4) 

When required, the function F(x) will be extended to the interval 
(-M) by 

F(-x) = F(x). (5) 

The function F(x) is, thus, the illumination required for a continuous 
aperture to produce the far-field pattern Hit). Equation (3) defines 
the array aperture as one, and the function Hit) to be bandlimited with 
bandwidth one. 

The approximation or synthesis problem consists in the determina- 
tion of the quantities X\ , • • • , x N subject to the condition of (1 ) so that 
S N (t) shall approximate H{t). 

In this form, the problem is that of numerical quadrature by means 
of an equal-coefficient rule. Sections II and III present methods for 
accomplishing this. Section IV drops the restriction of equal coefficients 
and applies the well-known Gaussian quadrature rule. Section V dis- 
cusses the existence of grating lobes and presents estimates for their 
location. 

II. A RIEMANNIAN SUM METHOD 

Let H(t) be a characteristic function, that is, H(t) satisfies the normal- 
ization condition (4) and the additional requirement 

F(x) ^ 0, < x < 1, (6) 

then the function 

L(x) = f F(u) du (7) 

Jo 

satisfies 

L(0) = 0, L(l) = 1 (8) 

and is monotonic increasing. Let 

y = L(x), x = G(y) (9) 

in which G iy) is the function inverse to L(x), then the required numbers 
Xj are given explicitly by 

Xi = G (cw 1 )' ls '' gff - (10) 
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The sum 

wo-^S-K 3 ^)] (11) 

is clearly a Riemannian sum for 

/ = f cos [tGiy)}dy, (12) 

and hence 

Km SM = I; (13) 

AT-»co 

however, 

/ = [ Fix) cos tadz = H{t) (14) 

and hence the approximation is secured. The error R N (t) given by 

R N (t) = H{t) - S N it) (15) 

will now be studied. For this purpose consider 

Lemma 1: ^ = f l gin [2imL(x)] sin tx dx 

Jo 

"♦ftr<0-4*£<-l) M T- 

ttN ik-i fc 

Proof: It will be convenient to introduce the function 

&r(fcf) = i g cos [«? (j, + ^WV\ > (16) 

thus 

OrftO) = 5,(0- (17) 

The function cos [tG(y)] may be expanded into a Fourier series on the 
interval (0,1); one has, 

cos [(Giy)] = Hit) + 2 a n cos 2tWJ/ + E &* sin 2imy, (18) 

n=l n=l 

in which 

On = 2 I cos [/(?(#)] cos 27r?iy dy, 

(19) 
b„ = 2 I cos [/(?(?/)] sin 2x71?/ d?/. 

•>0 



t 

a n = — c„ , 
xn 




, 1 — cos t 

b n = - 


1 A 
dn- 

irn 
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Define c„ , d n by 

c„ = / G'Uj) sin 2-imy sin [tG(y)] dy 

= / sin [2irnL(x)] sin tx dx, 

f , (20) 

d n = I G'(y) cos 2xny sin [tG(y)] dy 
J o 

= I cos [2tt7iL(x)] sin tx dx, 

then integration by parts applied to (19) yields 



(21) 



The Bernoullian function 

p(v) = \ - (v), (22) 

in which { y\ designates the fractional part of y, has the Fourier series 
,(,) = t *^». (23) 

n-1 7TW 

hence, replacing a„ , b„ in (18) by their values in (21), one obtains 

cos [tG(y)] = H(t) +LJT-cos 2vny 
t n=i n 

t " d (24) 

— - £ — sin 2irny + (1 — cos t)p(y). 
t „=i n 

By summation of the geometric series, one has 

^g««p[*-(»+V)]-^ ( -»""- *i* (25) 

= 0, iV | n, 

and hence, letting w = Nk (k > integral), 

~ g cos 2*-2\Tft (y + ^^) = ( " 1 )* cos 2t^%, (26) 
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| Xj sin MVl h/ + %^A = (-D* sin 2rtf%. (27) 

Equations (16), (24), (26), and (27) now yield 

S N (t,y) = Hit) +-LJT(-i) kC -^co S 2TNhj 
irN fc=i k 

- 4? S ( - 1 )* t* sin 27rJV% (28) 

7T7V a==i /c 

+ (l-cos,)iEp(.+ ^). 

The Fourier series for p(y), (23), permits ready establishment of the 
identity 

£4 + B w i ) = '( Ar » + 2)' (29) 



hence 



£„«,*/) = H(t) +-\ f T,(-l) kCj ^co S 2rNky 

TJ\ fc=l K 

- -4 I) (-D* nr sin 27riVfc 2/ (30) 

1TIV fc-1 AC 



Setting ?/ = in (30) yields the result of the lemma. 



. 1 — cos t 
+ —N- p 



Lemma 2: r ^ 2, integral, lF (r) (a) ^ £ r > or 

W M (x) ^ -£ r < for a ^ x ^ b 



f cos W ( 



aO d.r 



(r+l)/2„-(l/r) 



^ r2 irTI, "c; 



Proof: It is clear that only the inequality W (r) (x) ^ e r > need be 
considered. The case r = 2 will be considered first. The function W (x) 
is monotonic increasing, hence it vanishes at most once in [a,b], say at 
x = c, then 

/b re pb 

cos W(x) dx = J cos T7(i) dc + J cos IF(x) dx. (31) 
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Let ^ 8 S b — c be chosen, then 

•b r c+B 



f cos W(x) dx = [ cos W{x) dx + / cos W (x) dx, (32) 

Jc J c J c+S 



and hence 



One has 



r" I f b 

/ cos W(x) dx ^ 8 + / cos TF(x) dx 

Jc J c+S 



/ cos W(x) dx = I _-.. . d sin W(x-) 



+i TP'(*) 



= w „ ^T -, f dsinWix), 

W (C + 8) Jr+6 



W'(c + 5) Jc+s 

in which the second mean- value theorem was used, and hence 

2 



/ cos W(x) dx 



W'(c + 5) * 



Since 



-c+S 



W'(c + S) = [ W"(x) dx £ 8e 2 , 
one obtains, from (33), 

dx 



The choice 



yields 



/ cos W(x) 

8 = a/2 ef h 



^ s + u>- 



I cos W(x) 



dx 



^ 2 \/2 £2~ h . 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



The value of 8 in (38) may exceed b — c, however, in this case the 
inequality of (39) is certainly correct since the integral always admits 
the estimate b — c. 

Similarly choose ^ 5 ^ c — a, then 

/ cos W(x) dx = cos W(x) dx + / cos W(x) dx, (40) 

J a J a Jc-S 
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and hence 



f e I I f c ~ 6 

I cos W(x) dx\&\ I cos W(x) 

J a \ \ J a 



dx 



+ 8. 



One has 



rc-i rc-h j 

J cos W(x) dx = J ^7- 



a W'(X) 



d sin W(x) 



= wu- n I. dsinWr <*>. 



and hence 



Since 



,c-« 

/ cos W(x) 



dx 



< — 



W'(c - 8) ' 



-W'(c - 8) = f W"{x) dx ^ 5e 2 . 

Jc-6 

one obtains from (41 ) 



Hence 



and, from (31), 



/ cos W(x) dx 

I <* 

I cos W(x) dx 
cos W(x) dx 



*' + iS- 



^ 2 v^ £ 2 *, 



^ 4\^£ 2 i 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



The lemma is thus established for r = 2. 

Induction will now be employed. The lemma is assumed true for 
r = k ^ 2. Since TT ( * +1> (x) > 0, W lk) (x) is monotonic increasing, and 
hence vanishes at most once in [a,b], say at x = c. Choose ^ 8 ^ b — c, 
then 



[ cos TF(jf) dx = f cos W(z) dx + / cos lf(x) dx, (48) 
and hence 



/ cos W(x) dx ^ 8 4- / cos W(.r) da: 



(49) 
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The inductive hypothesis states 



hence 



Since 



one has 



The choice 



( cos W(x) di 

| Jc+i 

I cos W(x) 



g k2 w W w (c + 8)- im \ (50) 



^ 5 + k2 (k+1),2 W m (c + 5)- (1/fc) . (51) 



W m (c + 8) = f W (k+1) (x) dx ^ 8 £k+1 , 



(52) 



/ cos W(x) 



dx 



£ 8 + flS^'V^Vi • (53) 



. */2 -d/fc+1) 

= -S £fc+l 



yields 



/ cos W{x) 



dx 



^ (k + 1) 2 fc " £jt+1 



fc/2 „ -(1/fc+l) 



(54) 



(55) 



The inequahty of (55) remains correct even for 8 > b — c. 
Similarly, choose ^ 8 ^ c — a, then 

/•c rc—S re 

I cos TT(a;) dx = / cos JF(a;) d? + / cos W(x) dx, (56) 

Ja J a J c-s 



and hence 



dx 



I cos W(x) dx ^ / cos WXo?) rf: 

I •'a | I •'a 

The inductive hypothesis yields 



+ 8. (57) 



/ cos W(x) 



dx 



^ 8 + &2 (fc+1)/2 [-W (k) (c - 8)V llk) . (58) 



Since 



-W m {c - 8) = f W (k+1) (x) dx ^ 8e k+1 , (59) 
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one has 

I f cos W{x) dx I £ 5 + fc2 ( * +,) ' 2 r H/ *W (1/W . (60) 

I « 

Thus 

^ (& + l)2 W2 ft+ r (1/fc+l) , (61) 



/ cos W(.i;) tf.r 



and hence 



/ cos W(x) 

J n 



dx 



^ (k + l)2 {k+2>>2 £k+1 - (m+1) . (62) 



The lemma is now established. 
Theorem 1 provides an estimate of Rs(t). 

Theorem 1: r ^ 2, integral, L (r, (.r) ^ e r > or 

L (r) (.r) ^ -e r < /or ^ a; ^ 1 

=» | Ar(t) I £ r2 (r+1)/a - (1/r V- 1 - a/r) r(l + (l/r))er (1/r) 1 1 1 A r - 1 - (1/r> . 

Proof: One has, from Lemma 1, 

c„ = / sin [2icnL(x)] sin tx dx, (63) 

•'o 

and hence 

c n = \\ cos [2iroL(aO - tx] dx - i / cos [27rnL(.r) + te] <fc. (64) 
2 ^o 2 »'o 

Lemma 2 applied to the integrals of (64) yields 

|*|£ r2 <r+,)/2 - <1/r V- (l/r) £r - (1/r) n- (,/r) . (65) 

The infinite series for R N (t) in Lemma 1 may now be estimated. Using 
(65), one obtains 

| RAt) I ^ r 2 (H - 1)/2 - a/ 'V- 1 - (1/r V (1/r) I 1 1 N~ l - {,M E 71^177, . (66) 

Since the series of (66) is f (1 + (l/r)), the inequality of the theorem 
follows. 

An example of the above analysis is provided by the choice 

H(t) = J (t), (67) 
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that is the Bessel function of first kind and order zero. For this case 

2 
TrVl - X 2 ' 



^(*)=-^==5, (68) 



and hence 



Thus 



L(x) = -sm~ 1 x. (69) 



7T 



x = sin — y, (70) 



2 



and 



Xj = Sm 2-2N~- (71) 

The function L (x ) satisfies 

L'"{x) ^ 1 = *, (72) 

and hence, after numerical simplification, the error is estimated by 

\Rs{t) | < 8|*|i\T (4/3) . (73) 

Another example is given by 

H(t) = (smW¥) 2 - (74) 



One has 



and 



Thus 



Since 



F(x) - 2 - 2x, (75) 

L(x) = 2x - x 2 , (76) 

G(y) = 1 - VT^lj. (77) 

Xj = 1 - Vl - (2j - 1)/2A/. (78) 



L" = -2 = - £2 , (79) 

the error estimate obeys 

|JBjr(t) | < 1.3 U | N~ m \ (80) 
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If F (x) has high order of contact at the endpoints zero and one, then 
H(t) will decrease rapidly with increasing t, and hence the sidelobes 
will be small. In particular, let 

F U) (0) = 0, F u) (\) = 0, ^ j ^ k (81) 

then, integration by parts applied to (3) yields 

H(t) = -(-l) m W [ F lh+1) (x) sin tx dx, k even, (82) 

t Jo 



and 



H ( t ) = _(_!)<*-«/* 1 f F (k+1) (x) costxdx, kodd. (83) 

t Jo 

If F {lc+1) (x) is of bounded variation, then 

|F(f)|Sp, (84) 

in which V is the total variation of F ik+1) (x). Equation (84) shows the 
rapid decay of the sidelobes. 

An example of this type of tapered design is given by 

F k (x) = (2k + 1) (ff\ [rc(l - x)f (85) 

which has order of contact k — 1 and for which 

\H(t)\^V/t k . (86) 

In this case, V is the total variation of F k lk) (x). Since 

L« k+1) (x) = -(2k + Dl (*) = - £ *+* < ° < 87 > 

one has, from Theorem 1 , 

|Ar(f)| ^ EtltlN- 1 - 11 '*™, (88) 

in which E k is the constant determined by the theorem. 
The function 

H(t) = sm t/t (89) 

corresponds to 

F(x) = 1, L(x) = x, G(y) = y. (90) 
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The distribution of radiators is 

*/- (2i- l)/2tf (91) 

and therefore is uniform. The approximability of this function is poor 
compared to the previous examples. Theorem 1 does not cover this 
case since L" (x) = 0; however, the Fourier coefficients c n (20) may be 
explicitly evaluated, and the final determination of R»(t) obtained from 
Lemma 1. The result is 

«.(,) = « ±' f" ( ' + ff. , (92) 

■wN it=Z» k(t + 2irNk) 

in which the prime shows the absence of the term k = 0. Evaluation 
of the integral 



J P\ Nx + - J sin tx dx, 


(93) 


using (23), shows that 




Rtr(t) = — a? / P ( Nx + n) sm ix dx. 


(94) 


Since 




\p(Nx + |)| ^ 1, |sinte| ^ 1, 


(95) 


one has 




Iff* (01 ^ i\t\N~ l . 


(96) 



III. TCHEBYCHEFF QUADRATURE METHOD 

Let <p(x) e C (M) [— 1,1], then the Tchebycheff quadrature formula 1 is 

[ K(x)*{x) dx^l= Ev(«A ( K(x) dx = 1. (97) 
J-i M j=i J-i 

The fundamental points a } are determined by the conditions 

M [ x'K(x) dx = £ <*/ - fc , ^ ^ ^ ilf. (98) 

Define the polynomial a>(z) by the polynomial portion of the Laurent 
expansion of 



exp [Ml K(x) In (z - x) dx\ 



M 

= z 



( bi b 2 6 3 \ 



(99) 
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about the origin, 2 then the zeros of w(z) are the required numbers 
a,, • • • , a u . This procedure yields an approximation which, by (98), 
is exact if <p(x) is a polynomial of degree not exceeding M. To obtain 
an approximation to H(t) of the required form (1), one may set 

<p(x) = cos tx, K(x) = %F(x), M = 2N; (100) 

the points Xi, • • , x N are now chosen as those a, which are positive. 
Equations (3) and (97) yield the required result. 
Define the error, R N T , of Tchebycheff quadrature by 

R M T = f K(x)<p(x) dx-^Z *>(«,), (101) 

J-i M j-—i 

then Theorem 2 provides an estimate. 

Theorem 2: The real numbers «i , • • • , a M are determined as the zeroes 
of the polynomial co(rc) defined in (99) 

=>3-l<KnK./ = 4f K(x)o>(x)<p (U) (l;) dx. 

M ! J-i 

Proof: It will be shown that Tchebycheff quadrature is an instance of 
Newton-Cotes quadrature. 



Define 



>(*) 



[X — ajjco (a;) 
then the Lagrange interpolation formula is 

V(X) = £ *(«,) (M + *-^ «(*), (103) 

in which £ satisfies 

min (x, ai,---, a«) < £ < max (x, «i , • • • , ««)• (104) 
The coefficients of the Newton-Cotes quadrature formula are given by 

Cj = f K (x)(j(x) dx (105) 

and hence, one has 

[ K(xMx) dx=jt cMaj) + ± [_K(x)u(x)<p (M) (i;) dx. (106) 
Since Tchebycheff quadrature is exact when <p(x) is a polynomial of 
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Thus, secondary lobes may be produced of strength nearly equal to the 
main beam. These are called grating lobes. Their existence is the subject 
of Theorem 5. 

Theorem 5: There always exist grating lobes. 

Proof: Dirichlet's theorem* on simultaneous approximation states: 

Given Xi , • • • , x» , a positive integer q, and a positive integer r , 
there exists a number r in the range 

to^t^ r q N , (114) 

and integers p x , • • • , Vn , such that 

\rXi-Vi\ £ Vff. l^j^tf- (H5) 

Accordingly, choose r = 1 and t = 2itt, then 

txj = 2vtXj = 2irpi + (2x/g)0, | | £ 1, (116) 



and 



Thus 



cos toy = cos - > 1 - (27r 2 /g 2 ). (H7) 

9 



2x* 



&r(fl = £ Ay COS toy > 1 - -\ . (118) 

y=i 9 

Since g may be chosen arbitrarily large, the theorem is proved. 

An inspection of all the error formulae of this paper shows the common 
feature that they increase with increasing | t | and ultimately become 
trivial. For large 1 1 \, H(t) is small, hence, since S N (t), by Theorem 5, 
must ultimately become large, the error estimates must also become 
large. It follows that the grating lobe cannot occur until R N (t), RN T (t), 
or Rtr a (t) are at least one. The error estimates, therefore, provide a 
lower bound for the value of | t \ at which a grating lobe can occur. This 
is especially important in those designs where it is desired to eliminate 
the grating lobe from the scan sector. The methods of synthesis presented 
in this paper provide different estimates of location of the first grating 
lobe. Let T designate that location, then, for 

H(t) = J (t), T > $N il3 , (119) 

H{t) = (sin Wit)', T > .77N 312 , (120) 

H(t) = sin t/t, T > 2N. (121) 



- 



f 
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The above are the estimates obtained from the Riemannian sum method. 
The Tchebycheff and Gaussian quadrature methods do not yield esti- 
mates of grating lobe location nearly as advantageous as the Rieman- 
nian sum method. Thus, for 

H{t) = J (t), T > 4iV, 

(122) 
H(t) = sm t/t, T > 4jV. 

These results were obtained by rough approximations to the factorials 
in (112) and (113), however, they serve to show the difference between 
the Riemannian sum, and the Tchebycheff and Gaussian quadrature 
methods. Nonetheless, the last two methods may show a much smaller 
estimate of error for small I t I than the Riemannian sum method. 



REFERENCES 

1. Tomlinson Fort, Finite Differences, Oxford University Press, 1948. 

2. Milne-Thomson, L. M., The Calculus of Finite Differences, Macmillan and Co., 

1933. 

3. Sansone, G., Orthogonal Functions, Interscience, 1959. 

4. Titchmarsh, E. C, The Theory of the Riemann Z eta-Function, Oxford Univer- 

sity Press, 1951. 



